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The structure of the electron quantum size levels in spherical nanocrystals is studied in the 
framework of an eight-band effective mass model at zero and weak magnetic fields. The effect of 
the nanocrystal surface is modeled through the boundary condition imposed on the envelope wave 
function at the surface. We show that the spin-orbit splitting of the valence band leads to the 
surface-induced spin-orbit splitting of the excited conduction band states and to the additional 
surface-induced magnetic moment for electrons in bare nanocrystals. This additional magnetic mo- 
ment manifests itself in a nonzero surface contribution to the linear Zeeman splitting of all quantum 
size energy levels including the ground 15* electron state. The fitting of the size dependence of 
the ground state electron g factor in CdSe nanocrystals has allowed us to determine the appropri- 
ate surface parameter of the boundary conditions. The structure of the excited electron states is 
considered in the limits of weak and strong magnetic fields. 
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I. INTRODUCTION 

The modern technique of nanocrystal growth has 
created high outieal quality nanocrystals with narrow 
size distributiontia that can be also chemically doped.cl 
This allows one to study the level structure of|-jljiie| 
nanocrystals using intexband optical spectroscopyptl'El 
tunneling spectrosfiopy,0'El and the far-infrared intra- 
band spectroscopy.aLl Important additional information 
about the origin of these levels and their symmetry 
could be obtained by using optical magnetospectroscopy 
of the intraband and interband transitions and tun- 
nelling magnetospectroscopy, as was done many years 
ago in the case of atoms, hydrogenlike shallow defects 
and excitons. There is the other groucLjof experiments, 
such as magneto circular dichroism,tj spin-flip 
man scattering,Ell time dependent [Ea*aday rotation,E3't3 
EPR and ODMR jaeasurements,MOlla and magneto- 
photoluminescenceJlZI where the magnetic field is an es- 
sential component of the experimental technique. In- 
repretation of all these experiments requires the knowl- 
edge of electron and hole energy spectra in magnetic 
fields. 

A wide class of semiconductor materials can be pre- 
pared in nanocrystal form. In most cases they have a 
spherical shape and can be bare semiconductor nanocrys- 
tals or onion type heterostructure nanocrystals composed 
of an inner semiconductor core coated . w i th j«;veral spher- 
ical shells of different semiconductors Jl3ll3'E3 Technology 
allows one to vary the radius of the nanocrystal core 
from 10-40 A in a controlled manner. The energy spec- 



tra of the NCs are formed by the discrete quantum size 
levels (QSLs). The QSLs of the confined electrons in 
the conduction band are characterized by the value of 
the orbital momentum I and the total angular momen- 
tum I ± 1/2, similar to the bound electron states in 
atomsEJ. In a zero magnetic field, the levels are degener- 
ate with respect to the projection of the total momentum 
J- 

The spectra of the nanocrystals in an external mag- 
netic field H is well described as an atomic Zeeman ef- 
fect. A diamagnetic contribution to the spectra could 
be neglected because the nanocrystal radius, a, is consid- 
erably smaller than the magnetic length, L = ^ efi/cH, 
where e is the absolute value of a free electron charge, 
c is the speed of light. In a reasonable magnetic field 
H < IQT the diamagnetic contribution to the QSLs 
^ [a/L)^ ^ 1. However, the effect of the external 
magnetic field on the electron QSLs depends on the zero 
field spin-orbit splitting, Ac, of the levels with different 
j — I ±1/2 (One can find a description of the similar 
effect for atoms in Ref. |2l|). If Ac is larger than the Zee- 
man energy, the weak magnetic field splits the electron 
levels: 

AEf ^ fiBgfmH , m = -j, -j + 1, , (1) 

where — eh/2mQC is the Bohr magneton, mg is the 
free electron mass, m is the projection of the total mo- 
mentum on the magnetic filed direction and is the 
effective gr-factor of the corresponding state, that is an 
analog of the Lande factor for atoms. If Ac is small, 
the Zeeman splitting can be described as the sum of the 



2 



spin and orbital contributions: 



AEi = fiBQsSzH + fiBgdzH , 



(2) 



where — ±1/2 and Iz = —I, — / + I are the electron 
spin and angular momentum projections on a magnetic 
field direction correspondingly, and gs and gi are the spin 
and orbital g-factors of the electron in the corresponding 
state. In bare nanocrystals made of semiconductors with 
a simple parabolic conduction band gs is equal to the bulk 
electron effective g-factor gc and gi = nio/mc, where rric 
is the bulk electron effective mass. 

Spin-flip j-Raman scattering studies of CdS 
nanocrystalsjlil however, have demonstrated the 
dependence of the ground IS electron state g^ on the 
excitation energy, and thus on the energy of the electron 
state. The variation of gs with energy can be connected 
with iJae energy dependence of the bulk electron g- 
factor.Eil The eight band Kane model, that is used to 
describe the electron energy spectra in most direct gap 
semiconductors, gives the following de»eB4ence of the 
electron g-factor gc{E) on its energy, Et^ 



gc{E) = g* 



2E„ 



3 {Eg + E){Eg + A + E) 



(3) 



Here Eg is the band gap energy, A is the spin-orbit split- 
ting of the top of the valence band, Ep is the energetic 
Kane parameter, and g* = go + Ag, where go = 2 is the 
free electron g-factor and \Ag\ <C 1 is the contribution of 
remote bands to the bulk electron g factor. The second 
term in Eq. describes the negative contribution of 
the valence band (with A > 0) into the electron g-factor. 
It decreases with increasing the energy of QSLs caused 
by the reduction of the nanocrystals size. Thus, the en- 
ergy dependence of gdE) in Eq. (||) is conducive to the 
size dependence of gs in Eq. (H), that has been measured 
recently in bare and core-shell CdSe nanocrystals .113 

The size dependence of gs for the ground IS electron 
state in the spherical heterostructurc has been calculated 
in Ref. |2j within the eight band Kane model. It has 
been shown that in a spherical heterostructurc formed 
by two semiconductors, A and B, gs is the sum of the 
weighted volume contributio ns of e ac h mate rial and the 
interface contributions: gs = 5s (^) + gs{B) -f gAB- The 
interface term, gAB) is proportional to the square of the 
conduction band component of the wave function, /c, at 
the A/B heterointerface. This calculation shows clearly 
that gs is very sensitive to the value of the wave function 
at the heterointerface and to its leakage under the barrier 
and thus to the boundary conditions (BCs) imposed on 
the wave function at the heterointerface. In Ref. |2j 
the size dependence of the electron g-factor is calculated 
using standard BCs that assumes a continuity of /c at 
the interface. This leads to the vanishing of the interface 
term, gAB, in the bare semiconductor nanocrystals, that 
are modeled by the infinite potential barrier in the layer 
B because the standard BC at the surface in this case is 

/c = 0. 



However, the standard BCs are not always justified 
even for infinite potential barriers (see for example Refs. 

In general, the wave function /c does not vanish 
at the surface of bare semiconductor nanocrystals and 
satisfies the general boundary conditious. jLhat are the 
characteristics of the particular surf ace. Ej'lEI The effect 
of a semiconductor surface on the light absorption in in- 
direct band semiconductor nanocrystals has been studied 
in Ref. ^ and was found to be significant. 

In this paper we study the effect of the general bound- 
ary conditions on the QSLs of an electron in zero and 
weak external magnetic fields. We show that there 
is an additional surface-induced spin-orbit term for 
the electron states whose magnitude is proportional to 
A|/c(a)p . This surface-induced spin-orbit interaction 
leads the additional magnetic moment of the electrons 
in an external magnetic field in analogy with the addi- 
tional |£elativistic magnetic moment of the electrons in 
atoms .EZrEa This is conducive to the nonzero surface con- 
tribution to the electron g factor that controls the linear 
Zeeman splitting of the QSLs. The fitting of the exper- 
imental size |dependence of the electron g factor in the 
ground statdi3 has allowed us to determine the appro- 
priate parameter of the general boundary conditions for 
bare CdSe nanocrystals. The spectra of excited electron 
states has been considered in the two limits of the rela- 
tion between the zero field spin-orbit splitting and the 
electron energy in magnetic field similar to consideration 
conducted for atoms .til 

The paper is organized in the following way: In Sec. 
U we describe the eight -band Hamiltonians that include 
the external magnetic field effects. In Sec. IH we de- 



rive the general boundary conditions for the conduction 
band component of the envelope function at the abrupt 
surface (see Sec. [HA) and analyze the surface effect 



on the electron QSLs in a zero magnetic field (see Sec. 
IIIB). The effect of the weak external magnetic field on 
the electron QSLs is considered in Sec. [V, The size 



dependence of the electron ground state g factor in bare 
CdSe nanocrystals has been calculated and the value of 
the surface boundary parameter has b een determined by 
fitting experimental data in Sec. IVB. The symmetry of 
the electron excited states in an external magnetic field is 
studi ed in the low and strong magnetic field regimes (see 
Sec. [V C ) . The results are summarized and discussed in 
Sec. 0""^ 



II. EFFECT OF A WEAK MAGNETIC FIELD 
ON THE QSLS WITHIN THE EIGHT-BAND 
MODEL: PERTURBATION APPROACH 

A. The eight band model 

The energy band structure of the cubic semiconductors 
near the center of the first Brillouin zone caai. b e w ell 
described within the eight-band k • p modelBjfflEil In 
the homogenous semiconductor, the full wave function 
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ionS 



can be presented as the expansion 

Ai=±l/2 

E E K,(r)\Ra)u^, (4) 

where ui/2 and w_i/2 are the eigenfunctions of the spm 
operator ■S' — 1/2 ct, where <t = {^xi^yi^z} are the 
Pauh matrices, \S) is the Bloch function of the conduc- 
tion band edge at the F-point of BriUuoin zone repre- 
senting the eigenfunction of internal momentum / = 0, 
and = \X), \Ry) = \Y), \R,) = \Z) are the 

Bloch functions of the valence band edge at the F- 
point of Brillouin zone. The combination of these func- 
tions: 1/V2(|i?x) ± i\Ry)) and \Rz) are the eigenfunc- 
tions of the internal momentum I — 1 with projec- 
tions ±1 and on the z axis respectively (see Refs. 
|30| , ^ . The smooth functions ^'^^■^^2^'"') ^'^^ com- 
ponents of the conduction band spinor envelope func- 

tion = (3'^ y and n±i,2i'r). *^±i/2W> 



-1/2 



^z±i/2(''") ^^'^ ^iViZ components of the valence 



band spinor envelope vector 



*"l/2 
*"-l/2 



1/2 



J \ J \ ^2-1/2 / J 

ence of external magnetic field the eight-component en- 
velope function '^{r) = {^[r), \&^(r)} is the solution of 
the Schrodinger equationE3 



H{k) 



= E 



H''{k) ihPU2k 
-ihPU2k H"{k) 



(5) 



Here, the energy E is measured from the bottom of the 

conduction band, k = — (p+ - A) , where p = —iKV is 
n \ c J 

the momentum operator and A — (1/2) [H x r] is the vec- 
tor potential of the magnetic field. The 2x2 unit matrix, 
IJ2, and the Pauli matrices, Ox, <Jy and Oz, are acting on 
the spinor components of the wave functions (/x = ±1/2). 
P — —i{S\pz\Z) /rriQ is the Kane matrix element describ- 
ing the coupling of the conduction and valence bands 
(The Kane energy parameter: Ep = 2'mQP^, was intro- 
duced in Eq. (^, respectively). The off-diagonal matrix 
element in Eq. (||) acts on the as a scalar product 
- fe*" = kx'^Z + ky'i'l + kz'^z- The conduction band 
part of the Hamiltonian, H"^, acting on the spinor func- 
tion ^E"^ has the form: 



H^ik) 



ah 



1 



^^U2k^ + ^9*iiB{»H), (6) 



where g* is defined in Eq. (^, and a takes into ac- 
count the contribution of remote bands to the energy- 
dependent electron effective mass, mc{E). In cubic and 



zinc-blende semiconductors this dependence has the fol- 
lowing form: 



mc{E) 



E„ 



1 



En + E En 



E 



(7) 



In this paper we focus only on the electron QSLs with 
energies < E < Eg. This allows us to neglect hereafter 
the fc^ terms in the valence band part of the Hamilto- 
nian, HZj_and to present it in the spherical approxima- 
tion asMM 



H^ik) 



Eg + -A + {1 + 3K) {IH) U2 



^/\{ia) + ]^^iBgo{^H). 



(8) 



Here, the Hamiltonian H"" should be considered as the 
2x2 matrix acting on the spinor vector rather 
than the 6x6 matrix Hamiltonian acting on the six 
component wave function as in Rcf. ^ Correspond- 
ingly, / = {ix,Iy,iz} in Eq. (||) is the vector oper- 
ator. It is straightforward to show that in this repre- 
sentation (/T)*" ^ -i[T X where T is an arbi- 
trary vector. The terms proportional to A in Hamil- 
tonian m) describe the-effect of the spin-orbit coupling 



Epl&EgfB where is the magnetic Lut- 



^ One can find that k = —2/3 from 
= (-2 - 7f + 271- +-373'^)/3, where 
the effective mass Luttinger parametersj23 j^^2 3) can be 
presented as: 7^ « Ep/SEg atui Ep/QEg in 



and K = K — 
tinger paraip 
the relatione; 



the approximation used above .CJ 



B. Linear Zeeman effect as a perturbation 

Following the perturbation approach developed in Ref. 
|, we present the Hamiltonian (||) as: 



HpBik)^H%p) + H {p,H) , 



(9) 



where the Hamiltonian H'^ describes the electron states 
in zero magnetic field: 

= ( ^""^P" - ^^^^^ . ^ (10) 

1^ -^PU2P iA (la) {Eg + iA) U2 I ' ^ ^ 

and the Hamiltonians H describes the effect of a weak 
external magnetic field. Linear in the magnetic field 
terms of H (p, H) can be written: 

H {p, H) = hbHG + -{A, V} = hbH{G + G), (11) 
c 

where {a,b} — l/2(a6 + ba) denotes an anticommutator 
and V = dH^{p) j dp is the envelope velocity operator. 
The operators G and G are defined as: 



^g*ian) 
' iin)U2 + yoi^n) 



(12) 
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G 



a(£n) 



mo , 



■ r> ^ r 
i_r-— [ 







U2 



(13) 



where n = H/H is the unit vector directed along the 
external magnetic field, and i = |- [r x p] is the envelope 
angular momentum operator. 

Generally, the calculation of the linear Zeeman effect 
for the carriers confined in a 31? external potential of 
arbitrary shape is a straightforward theoretical problem. 
The same procedure is used for an impurity potential and 
for a zero dimensional quantum dot of an arbitrary shape. 
Using the unperturbed Hamiltonian one has to find 
first the energy levels E\ and the corresponding eigen 
functions \1/a = {^'^j of the carrier in the external 
confining potential. Then, the effect of the magnetic field 
can be found perturbatively: 



(14) 



where 
obtain: 



the wave functions should be normalized 
i = 1. Substituting ^'a into this expression we 



A^A = AEc + AEy + AE,oup , (15) 

where the contribution of the valence and conduction 
bands to Zeeman effect AEc and AE^ correspondingly 
are: 



(*1.go(^n)|*") + [^l\{In)\^l 
and the coupling correction is given by 

AEcoup / ,T,c I ■ "^O-P r II ,T. u 

= i^h — - — C/2 n X r * 



(16) 
(17) 



(18) 



In principle, this straightforward approach allows us to 
find the linear Zeeman splitting of the energy levels in 
any arbitrarily shaped hetero-nanostructure. However, 
the calculations of the zero field wave functions ^\ for 
arbitrary shape heterostructure is a rather cumbersome 
procedure and, consequently, the Zeeman splitting cal- 
culation that uses these wave functions is very compli- 
cated. In the present paper, from hereafter, only spheri- 
cal semiconductor nanocrystals are considered. The high 
symmetry of these quantum dots allow us to calculate 
the Zeeman splitting of the QSLs. 

III. ENERGY LEVELS OF THE ELECTRON 
CONFINED IN BARE SPHERICAL 
NANOCRYSTALS 



General boundary conditions for the envelope 
wave functions 



We will consider bare semiconductor nanocrystals 
which surface can be modeled by an impenetrable bar- 



rier. To find unperturbed envelope wave functions that 
are described by the Hamiltonians of Eq. ( 10 ) in the bulk 
region one should impose an appropriate boundary con- 
ditions on these functions at the nanocrystal surface. For 
the impenetrable barrier the general boundary conditions 
(GBCs) have been shown in Ref. ^ to guarantee the 
vanishing of the normal to the surface component of the 
envelope fiux density matrix. In a spherical nanocrystal 
the normal envelope flux density matrix, J^^{r), should 
vanish at any point of the nanocrystals surface: 



A/3 I 



T • y*A,*, 



0, (19) 



where a is the nanocrystal radius and t — r jr (in spher- 
ical coordinates r is defined as r = (r, 0,(p)). This 
general requirement of Eq. (|l^) should be satisfied for 
two arbitrarily chosen eigenfunctions ^\ and of the 
Hamiltonian defined in Eq. (|l^) with energies E\ 
and respectively. Substituting the explicit expres- 
sions for the normal components of the envelope veloci- 
ties T-V = T- dH^/dp into Eq. (|l|) one obtains: 



2moP 



ih 

2mn 



[a {%T ■ VVI/5 



h 



(20) 



Using equation = E^ one can express the compo- 
nents of the spin vector SP" through V^P^: 



+ 



2moP 

ih 
AmoP 



mo 



mdE)^ 
iUE) - g*) [axV , 



(21) 



where the energy dependent electron effective mass 
mc{E) and g-factor gc{E) are defined in Eqs. (0) and 
(^) respectively. The normal projection (x^") can then 
be written: 



T*''(r) = 



2moP 

h 

■inioPr 



mo 



mc{E) ) dr 
Ag{E){aL)-^^ 



(22) 



where Ag{E) = g* - gc(E). Substituting Eq. (g2|) into 
Eq. ( pO| ) one obtains the general requirement for the con- 
duction band component of the envelope wave function 
at the NC surface: 



5 



ZniQ 



^0 ,..cdn 



mo 



dr 



= 0. 



(23) 



Equation (^) must be fulfilled at each point of the sur- 
face. It is satisfied for all arbitrary chosen energy states 
A and 77 if and only if the boundary conditions for as 
well as for 'J'^j are given by: 



^'ir = a,e,^) = Tao 



Too (95''=(r, e,(/7) 



dr 



2r 



Ag{E)ia-L)^''{r,Q,ip) 



(24) 



where T is a real number constant independent of E and 
flo is the lattice constant. We assume that the surface 
of the nanocrystal also possesses the spherical symmetry 
and is characterized by the same BCs at any point of the 
surface. In this case parameter T in the GBCs of Eq. 
(U) is independent of angles 9 and 



B. Electron energy levels: the surface induced 
spin—orbit splitting 

Let us consider the effect of the GBCs of Eq. (p^ ) 
on the electron QSLs in the absence of a magnetic field. 
Using Eq. ( pT| ) that expresses the valence band com- 
ponent of the wave function through its conduction 
band component ^''^ one can derive the bulk Schrodinger 
equation describing inside a semiconductor nanocrys- 
tal: 



2TOe(-B) 



C/2V^*^(r) = E^ii^ir) . 



(25) 



One can see that each spin component of ^f"^ is a so- 
lution of the standard bulk Schrodinger equation with 
the energy dependent effective mass mc{E). Equation 
(^) does not contain a spin-orbit term and its solution 
can always be presented as fi{r)Yi^i^{Q, tp), where I and 
I2 are the angular momentum and angular momentum 
projection, correspondingly, Yi^i^ {Q, if) are the spherical 
harmonics, and the radial function fi{r) satisfies to the 
following equation: 



A, 



2meiE) 

92 2d 



Aifi{r)^Efi{r), 
1(1 + 1) 



(26) 



Q2j. J, Qj, 



At the same time the GBCs of Eq. ( p4| ) mixes two 
electron spin sublevels, and all electron states are charac- 
terized now by the total angular momentum: j — L + S. 



As a result is an eigen spinor of the operator of the 
total angular momentum: 



vI/'= = /±(r)%,,™(e,vp), 



(27) 



where f2j.;^m(0, <^) are the spherical spinors, j and m are 
the total angular momentum and its projection respec- 
tively, and I is the angular momentum of the electron 
state. The two radial functions fj^ describe the two elec- 
tron states with j = lzLl/2 (note that the states described 
by /;~ exist only for / > 1). We use the definition of the 
spherical spinors as in Ref. Isl 



'-'l/2,l/2;Z,m-l/2-^',™-l/2 
'^l/2,~l/2;;,m+l/2-^',™+l/2 



(28) 



where j = I ±\/2 and Cji"^i.j2 m2 '^^^ Clebsch- 
Gordan coefficients. Substituting Eq. ( p7| ) into Eq. ( |2^ ) 
one can write the general boundary condition for the ra- 
dial wave functions (r) : 



/±(a) = Tao 



Too 
nirXE) 



f^'{a)+5kfff{a) 



where 



I 

2a' 



Sk+^—AgiE), ; = 0,1,2..., 



for the states with j = I + 1/2 and 



6k^ = 



l + l 
2a 



Ag{E), Z = 1,2,3. 



(29) 



(30) 



(31) 



for the states with j = 1 — 1/2. The conventionally used 
standard BCs for the impenetrable barrier assumes van- 
ishing of the wave function at the nanocrystal surface, 
and thus correspond to Tag = 0. In general, however, the 
conduction band component of the wave function does 
not vanish at the surface. Furthermore, it has been shown 
in Ref.^ that the boundary condition with Tag = is 
incorrect in cases where inter band coupling is important. 
The appropriate value of Tao should be determined from 
the fitting to the experimental data. 

The solutions of Eq. ( p6| ) can be written as f^{r) — 
[C^ / a'^/'^)ji{(j)f^r / a), where ji are spherical Bessel func- 
tions, Cf^ is the normalization constant determined by 



the condition /(|*''P + 
nected with the energy. 



ib'"\^)dV = 1, and (t)f^^ is con- 
E, of the electron level: E = 
h'^(f)f^/2mc{E)a'^ . An equation that determined (t>fn 



6 



obtained by substitution of fj^{r) into the BCs of Eq. 



4mn 



Ag{E)^" {aL)^'' ] -dS , (33) 



(32) 



Tap mp : 



I 



21 + 1 



l + i 

21 + 1 



Ln 



The n solution of this equation defines the energy of the 



,th 



electron level with the angular momentum I and the normalization condition: 



where dS = r^sin{Q)dQd(p. Note that although '^'^ 

= the normalization constants Cf ^ Cf . 

However, in the considered approximation they can be 
replaced with high accuracy by the averaged constant 
Ci = C;"*" = C;~ that is determined by the approximate 



/ |/i|Vdr+ / |$j|2r2dr = l 
Jo Jo 



(34) 



total angular momentum j = I ± 1/2. In the case of the 
standard BCs Tap — 0, and the solution is the n-th 
zero (pf „ of the spherical Bcssel functions ji . The general 
BCs with Tap ^ takes into account the effect of the 
surface on the electron energy levels that is important = 
in small nanocrystals. In large nanocrystals satisfying 
the condition a 3> \Tao\(rnQ/mc), the solutions (j>^^ are 

close to (h° and the effect of the surface on the confined t^^at neglects the second term in Eq. (|l|). Thus, if 



mp \ /, /|2 , '(^ + I)|^j2 
2TnQEp V" mc{E) ' ^ '''' 



■\fi\ 



electron states is negligible. 

The electron states with j = I + 1/2 and j = I — 1/2 
have different energies for / > 1 as a result of the general 
BCs of Eq. (29). This difference describes the surface 
induced spin-orbit splitting of the excited states: Ac = 
Ei+i/2 — Ei~i/2- It can be shown that Ac is positive if 
the spin-orbit splitting of the valence band A > 0. 

Figure [l] shows the effect of the general BCs on the 
energy of the ground IS' {I = 0) and the first excited IP 
{I = 1) electron states in bare CdSe nanocrystals. The 
surface induced spin-orbit interaction splits the last state 
into two states with the total angular momentum j = 3/2 
(P3/2 state) and j — 1/2 {P1/2 state) and the size depen- 
dence of this splitting Ac = i?(lP3/2) — £'(1^1/2) is shown 
in Fig. |l|(b). The calculations have been performed for 



the following bulk parameters of CdSe: Ep = 19.0 eV 
and TOc(O) = 0.116 mp from Ref. H, Eg ^ 1.839 eV, 
A = 0.42 meV, and gc(0) = 0.68 from Ref. ||, that 
result in a = —1.07 and g* — 1.96. One can see that 
even a small surface parameter introduced by the GBCs, 
I Tap I = 0.6 A, significantly affects both ground and ex- 
cited states in small nanocrystals. The positive (nega- 
tive) value of the surface parameter Tap > (Tap < 0) 
increases (decreases) the energy of all states in compari- 
son with Tap — 0. The relative order of the P3/2 and P1/2 
states with the total momentum j = 3/2 and j = 1/2 
always coincides with the relative order of the valence 
band subbands characterized by the same total momen- 
tum. The value of the spin-orbit splitting, Ac, however, 
remains for all excited states much smaller than the av- 
eraged energy Ei = [lEi^^i^ + (l + l)Ei+i/2) / {21 + 1) of 
the level with the orbital momentum I. 

If \6kfTaQ\ <^ 1, the averaged energies Ei and corre- 
sponding wave numbers can be found from the sim- 



plified BCs given by Eqs. 



3) with Skf 



0. 



The 

can be found perturbatively. In the Appendix we obtain 
the following energy corrections: 



energy corrections AE^ coming from the small Skf^ 



AET = 



2mp 



\Sk^Tao\ ^ 1, Ac < El, the energy levels of the ex- 
cited states Ef with j = 1 + 1/2 are the eigen energies 
of the effective Hamiltonian Hi : 



Hl^l±l/2d,m 

Hi = E1U2 + 



El ^l±l/2. 



I An ■ 



1 



21 + 1 



A,{Ei)[La 



(35) 



where the energy dependent spin-orbit splitting Ac(£'i) 
is given by 



A,{Ei 



1 + 



1 



2 J 2rnpa- 



:^sur{El), (36) 



^sur{Ei) = Ag{Ei)a''\Ji{af = A.g(£;OCf jf (0;,„) . (37) 



The approximate expressions of Eqs. (36,3^) describe 
the exact splitting Ac(-Ei) of the IP states in CdSe 
nanocrystals shown in Fig. |^(b) with high accuracy (the 
largest error in the smallest nanocrystals is about 2% for 
Tap — —0.6 A). It is important to note here that while 
the matrix elements of the spin-orbit operator {La) are 
nonzero only for Z > 1, Eq. ( |37| ) defines also the surface 
parameter Yigur for the I = Q S symmetry states. 

The dimensionless parameter Yigur is proportional to 
the spin-orbit splitting in the valence band A and to the 
square of the wave function at the surface. It describes 
the surface-induced spin-orbit coupling of the conduc- 
tion electron QSLs in bare semiconductor nanocrystals. 
In bulk semiconductors the integration in Eq. (p3) re- 
sults in zero spin-orbit splitting of the conduction band 
states because it must be carried out at the remote 
bounding surface where the wave function is vanishing. 
Ssur = if one assumes the vanishing of the conduction 
band component wave function at the surface of spher- 
ical nanocrystals. This assumption is never justified if 
the coupling between conduction band and valence band 
components is significant. Thus, the spin-orbit splitting 
of the electron QSLs is caused by the admixture of the 
valence band states near the surface. 
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IV. THE FINE STRUCTURE OF ELECTRON 
QSLS IN A MAGNETIC FIELD 

A. Surface induced magnetic moment of the 
confined electrons 

Let us consider now the energy of the electron QSLs, 
^ in a weak magnetic field. Substituting the zero field 



wave functions = {5''^, '4'"} given by Eq. (27) and 
( pf] ) into Eqs. (|l^) and (|l^), we obtain for the energy 
correction AE fvAEr + AErouv- 



AE- 



1.1 bH 
mo 



1 



^^ME)ian)\^f^) + 



(in)|*" 



(38) 



where we neglect AEy corrections of Eq. ( |17| ) as being 
smaller by the factor Ei/{Ei + Eg). The second term of 
Eq. ( p^ ) describes the unexpected surface contribution 
to the electron magnetic energy: 



H 



J dS^Ag{Ei){'^''*[r x [a x r]]-^') .(39) 



The effective Hamiltonian describing the fine structure 
of the electron state with the orbital momentum I in a 
magnetic field can be written as Hi + Hh, where Hi is 
defined by Eq. (|3^) and the effect of a weak magnetic 
field is described by: 

Hh - ^^lB9s{El){&H) + fiB9i{Ei)iLH) + 

{El) {[t x[&x t]H) . (40) 

Here the weighted spin Qs and orbital gi g factors are 
given by 



gs = / gc{Ei)r'dr\fi{r)\\ 



(41) 



gi = I mo/m,iEi)r'dr\fi{r)\' , (42) 

'0 



respectively. The first two terms in Eq. ( [fOD describe the 
averaged volume energy of the spin ^Jig — —^Bgs{Ei)S 
and orbital /x^, = ^lt-B9i{Ei)L magnetic moments in 
an external magnetic field, and the last term describes 
the energy of the surface-induced magnetic moment 
(J'suriEi) given by: 

t^sur = - ^ ^sur {El )[t X [& X t] . (43) 

This magnetic moment arises from the surface-induced 
spin-orbital term of Eq. (35) which is modified by the 
Larmor precession of the electron in an external mag- 
netic field. Indeed, the surface energy term AEsur — 
— {fig.^j,H) can be obtained directly by replacing — zV — > 
— iV + {e/ch)A in the spin-orbit perturbation term 



{La) = -i{[r X V]d-) of Eq. (|3^). The origin of the 
surface-induced electron magnetic moment is similar to 
those of an additional relativistic magnetic moment of 
the electron in atoms, and their values are derived in a 
similar way (see for example Rcfs. |2T| , p7| , p8| ) . It is neces- 
sary to note that the surface-induced magnetic moment 
fj-g.y^r does not vanish for the states with I — even when 
it arises from the spin-orbit coupling term {La). 



B. Size dependence of the ground state electron g 
factor 

All electron states that have S symmetry are Kramers 
doublets that are degenerate with respect to its spin pro- 
jection. The external magnetic field lifts this degeneracy 
and splits these states into two levels with energy: 



Eo{H)^Eo±^^iBgs{Eo)H, 



(44) 



where the ± signs correspond to the electron states with 
spins parallel and antiparallel to the magnetic field di- 
rection. The spin effective g factor can be obtained from 
Eq. M) as: 



gs{Eo) = gs{Eo) + 



(45) 



The effective electron g factor for the S electron states 
in spherical heterostructures, gs{Eo), was first obtained 
in Ref. ^ where the standard BCs at the heterointer- 
face were used to describe the size dependence of the 
electron g factor. In the case of the bare semiconductor 
nanocrystals the standard BC at the surface {fo{a) = 0) 
would lead to zero surface contribution to the effective 
electron g factor. We have examined here the effect of 
the general BC of Eq. ( p9| ) on the electron effective g 
factor in CdSe nanocrystals. The weighted g factor of 
the states with S symmetry in Eq. (B5h can be written: 



UE)^9o + / r^dr\fi{r)\^g,{E)-go) 



(46) 



where the corrections AEy of Eq. ( p7| ) to the (js of Eq. 
( PI ) have been added. Figure g shows the size depen- 
dence of the ground state electron g factor, gs{Eo), cal- 
culated using the standard BC with \Tao\ — and the 
general BC with \Tao\ = 0.6 A. One can see that the ex- 
perimental size dependence of the electron ji-factor taken 
from Ref. |l^ is described very well by the general BC 
with the negative parameter Tap = —0.6, while the use 
of the positive parameter brings the theoretical curve far 
away from the data. 

The size dependence of the bulk-like energy- 
dependent g factor, gc{Eo), calculated with the help of 
Eq. (H) for Tao = is also shown in Fig. ^. One can 
see that the bulk g-factor gc fails to describe the ex- 
perimental data even in the largest dots. At the same 
time, the difference between two gs curves calculated 
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with Tqq = —0.6 A and Tag = is not very large. 
This is because the general BCs affects the electron g- 
factor in two ways: indirectly, through the change of 
the zero field energy Eq decreasing/increasing the elec- 
tron g-factor, and directly, through the surface contribu- 
tion l/STisuri^S), that is always positive. In the case 
of the negative Tqq < the GBCs decreases the energy 
of the ground electron state decreasing the electron g- 
factor, and these two effects partly compensate for each 
other. However, the surface contribution l/3Es„r(l'S') 
calculated with Tap = —0.6 A is significant in small 
nanocrystals (see Fig. ^). 



angular momentum I. Assuming that the states with dif- 
ferent / are not mixed, one can find: 



nr. 

I2 



{2l + 3){2l~l 



{2l + 3){2l-l) ""^ 
■{L,Lj}]Yi^m{Q,^), (51) 



where i, j — x,y, z and dij is the Kronecker delta symbol. 
Substituting Eq. ^ into {H[t x [ct x x]) = {Ha) - 
{HT){Ta-), we rewrite the magnetic field Hamiltonian 
Hh of Eq. (|o|) as: 



C. Symmetry of the confined electron excited 
states in an external magnetic field 

The fine structure of the excited electron states in an 
external magnetic field is very sensitive to the relation 
between Ac{Ei) and the energy of the orbital magnetic 
moment gifisH. In the low-field regime, when the mag- 
netic energy is smaller than the separation between the 
electron states with j = I zL 1/2, the external magnetic 
field satisfying the condition giusH <C Ac splits the lev- 
els according to projection m of the full momentum j on 
the magnetic field as: 



E 



Ei+Ac{Ei 



I 



21 + 1 



tiBgt{Ei)mH, (47) 



for j = / + 1 /2 and 



^Ei- A,{Ei)^±^ + MS g;{Ei)mH. (48) 

for j = / — 1/2. The effective electron ^-factors 

2j + l 



+ 1 .X 2i - 1 , , 



2.f 



2j 



8J(J + 1)^^"'"' ^^^^ 



9, 



gs[Ei) + ———gi[Ei) - „ . , _,, Ss«r 



2(J + 1) 



2(J + 1) 



8j(j + 1) 



(50) 



are the analogs of the Lande factors for electrops in atoms 
in the case of the "anomalous" Zeeman effect .til The sur- 
face contributions to the electron (/-factors (oc S^ur) are 
similar to the relativistic corrections in Refs. |27| , p8| . 

When the zero field fine structure splitting becomes 
smaller than the energy of the orbital magnetic moment, 
the magnetic field Zeeman term mixes the states with 
different j = / ± 1/2. In this case the projections of the 
spin momentum Sz and orbital momentum Iz on the di- 
rection of the magnetic field are a more convenient nota- 
tion for describing the electron state fine structure. The 
spin-orbit term (LS) and the additional surface moment 
(figy^^H) may, however, mix the states with different val- 
ues of the Szlz product. 

Let us consider the matrix elements of the surface mag- 
netic moment operator /Xju^ on the eigen- functions of one 



Hh = ^iB9s{El){SH) + fiBgi{Ei){LH) + 
,^^ {SH)L^ + {{SL){LH)} 
^^^^"'■^^'^ (2Z + 3)(2Z-1) 

Here the spin electron g factor is given by 

1 



(52) 



gs{Ei)^UEi)- 



{2l + i){2l~l) 



(El) . (53) 



The matrix elements of the last term in Eq. (|5^) are 
zero for / = 0. If I^suriEi) <^ gi{Ei) one can neglect the 
off diagonal elements of the operator {{SL){LHy\ and 
replace them with {SH){Ln)'^. This is the case for the 
first excited IP state in CdSe nanocrystals (for compar- 
ison, see the corresponding curves in Figs. ^ and ^). In 
this case the last term of Eq. ( |5^ ) describes the surface 
contribution to the spin splitting of the electron levels 
that depends on the pmjection of the orbital momentum 
on the magnetic field. EZl 

In the strong field regime (similar to the case of the 
"quasi-normal" Zeeman effect or "copiplete" Paschen- 
Back effect for the electrons in atomsO) one can addi- 
tionally neglect the off diagonal elements of the spin- 
orbit operator (LS). The fine structure of the electron 
level with orbital momentum I is described by: 



Ei,i^,s. =Ei + ^^A,{Ei) + 

IJ-B Szgs{ElJz) H + lJLBlzgi{El) H , 



(54) 



where Iz ~ —l,...0,...l and Sz = ±1/2 are the projec- 
tion of the electron angular momentum and spin on the 
magnetic field direction and 



gs{EiJz)^gs{Ei) 



p+i+n 



1 



(2/ + 3)(2Z- 1) 



\Ei). (55) 



Thus the structures of the excited QSLs with / > 1 
may be different in weak and strong magnetic fields and 
undertake the anticrossing in the intermediate magnetic 
field (similar to those known for atom£d). In nanocrys- 
tals the weak, intermediate and strong field regimes de- 
pend strongly on the nanocrystal size as well as on its 
surface conditions. 
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V. DISCUSSION AND CONCLUSION 

We have studied the effect of general boundary con- 
ditions (the surface effect) on the electron QSLs and 
their Zeeman splitting in spherical bare semiconductor 
nanocrystals. The above consideration that has been 
carried on in the eight-band effective mass model can 
be easily extended to describe spherically layered het- 
erostructures and wide gap semiconductor nanocrystals, 
which are better described by the fourteen-band model 
(see for example Ref. |3^ ). 

Comparing the results of our theoretical calculation 
with the experimental size dependence of the electron g- 
factor we have determined the surface parameter Tag — 
—0.6 ±0.05 A in bare CdSe nanocrystals. Fitting the ex- 
perimental data, we used Ep = 19.0 eV, that has been ob- 
tained independently from bulk measurementsE3 and that 
describes better the experimental data than Ep — 17.5 
eV, previously used for CdSe in Refs. The surface 

parameter Ta^ characterizes the electronic properties of 
the surface and should be considered as the additional 
one to the set of parameters that describes the bulk prop- 
erties of semiconductors. In CdSe nanocrystals prepared 
by a different technique Tap can be different. 

The extracted absolute value and the negative sign 
of Too = —0.6 A is consistent with our theoretical ex- 
pectations for Tag in studied CdSe nanocrystals. Its 
value is very close to the theoretical value of the sur- 
face parameter for semiconductors with a symmetrical 

band structureEl \{Ta^)s\ = a* = ■<Jh^/2Epmo ~ 0.45 

A. It can also be shown within the eight band effective 
mass model, that the negative sign of the surface param- 
eter does not allow the existence of the surface localized 
states with E < (gap states) (one can find similar con- 
sideration in Ref. E6h. Indeed, the CdSe samples studied 
m Refs. [l||l| show very high PL quantum efficiency 
and do not show deep gap transitions. The negative pa- 
rameter Tao leads also to an additional nonparabolicity, 
bowing the size dependence of the electron energy levels 
(see Fig. ^ and therefore may describe the unexplained 
experimental size dependence of the IS electron level in 
small CdSe nanocrystalsE 

Using the GBCs we have found also the direct surface 
contribution to the spin-orbit effects in zero and weak 
external magnetic fields. The surface contribution to the 
zero field spin-orbit splitting is similar to the interface 
contribution obtained in Refs. |3^,^ for 2D electrons 
in planar quantum wells by using the spin-dependent 
boundary conditions. It has been pointed out in Ref. ^ 
that the interface contributicm-to the Rashba spin-orbit 
term in the 2D Hamiltoniancil is related to discontinu- 
ity of the band parameters at the semiconductor het- 
erointerface and that this contribution is an additional 
one to those connected with the space charge and/or the 
external electric field. The importance of the effects de- 
scribed by the Rashba tjeijipa.for the 2D electrons confined 
near the curved surfaceBa and cylindrical semiconduc- 



tor quantum dotsB has been emphasized recently. The 
Rashba spin-orbit term in our spherical dots is a direct 
consequence of the GBCs for the envelope function. The 
same consideration can be made for cylindrical dots or 
any other nanostructure geometry (the results will be 
published elsewhere). The GBCs provide an important 
connection between the constant that describes the mag- 
nitude of the spin-orbit term and surface conditions in 
nanostructures. 

In conclusion, we have shown an important influence of 
the semiconductor surface on the electron energy struc- 
ture in bare spherical nanocrystals. The effect of the sur- 
face has been modeled through the choice of the bound- 
ary condition parameter that describes the nonzero value 
of the envelope function at the nanocrystal surface. The 
additional nonparabolicity of the quantum size energy 
levels, the spin-orbit splitting of the electron quantum 
size levels, and the additional magnetic moment of the 
electrons have been shown to be induced by the sur- 
face. The effects are significant in small nanocrystals and 
their considerations require a multiband effective mass 
approach because interband coupling is important there. 
The analysis of the experimental data allows us to deter- 
mine the appropriate parameter of the boundary condi- 
tions that characterize the surface in studied bare CdSe 
semiconductor nanocrystals. 
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APPENDIX: EFFECT OF THE SMALL SURFACE 
PERTURBATION ON THE QUANTUM SIZE 
LEVEL ENERGY IN SPHERICAL 
NANOCRYSTALS 

We are interested in deriving a variation of the electron 
level energy E caused by a small variation of the parame- 
ter <j that characterizes the surface BC for the bulk wave 
function /(r, <j). The GBCs at the spherical surface of 
the nanocrystal with radius r — a can be written as 



/ (a,0 = /(«,0^(^), 



(A.l) 



where A(<r) is the real number constant. The func- 
tion /(r, <;) is the solution of bulk Schrodinger equation 
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Dk'^f{r,<,) = E{(,)f{r,<,) where the constant D is in- 
dependent of q. Taking a derivative of the Schrodinger 
equation on multiplying both parts of the resulting 
equation by /* and integrating it over the sphere vol- 
ume, one can obtain: 



ergy variation: 



(A.3) 



dE 



= Da-" 



df 



df 



and Ais) = 



Taking a derivative of Eq. (A.l) and substituting rection AE^ 
df'ja <;)/d<; = dA{^)/d^f{a,<i) + A{<;)df{a,i)/dc; into 
Eq. (]A.2|) one arrives at the final expression for the en- 



(A.2) Substituting D 

mc{Ei)/mQ[l/ (Tao) + c;] with ^ = -~Skf{Ei) into 
this equation we obtain Eq. (33) for the energy cor- 
to the averaged energy of the electron 



level with the orbital momentum I caused by the small 
perturbation of the boundary conditions of Eq. (E9) . 
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FIG. 1: The effect of general boundary conditions on the ground 15 and first excited IP energy levels (a), and on the spin- 
orbit splitting between IP3/2 and IP1/2 levels Ac (b) in bare CdSe nanocrystals. The size dependencies are calculated for the 
standard BCs with the surface parameter Tao = A (solid line), for Tao = —0.6 A (dashed lines), and for Too = 0.6 A (dotted 
lines). The bulk parameters of CdSe used in calculations are described in the text. 



FIG. 2: The size dependence of the ground IS state electron g- factor in bare CdSe nanocrystals calculated for the standard 
BCs with the surface parameter Tao = (solid line), for Tan = —0.6 A (dashed line), and for Tao — 0.6 A (dotted line). 
The empty squares show the experimental data from Ref. uA. For comparison we also show the energy dependence of the 
bulk electron g-factor gdE) calculated at the energy of the IS electron level (dash-dotted line) and bulk value of the electron 
g-factor (short dash-dotted line). 



FIG. 3: The size dependence of the dimensionless surface parameter Esur of Eq. ( pTI) for the ground 15 (/ — 0) and first 
excited IP {I = 1) states in bare CdSe nanocrystals. The dashed and dotted curves are calculated with the surface parameters 
Tao = —0.6 A and Tao = 0.6 A respectively. 



FIG. 4: The size dependence of the orbital g- factor, gi {lP), for the first excited state in bare CdSe nanocrystals calculated for 
the standard BCs with surface parameter Tao — (solid lino), for Tao — —0.6 A (dashed linos), and for Tao — 0.6 A (dotted 
lines). For comparison we also show the bulk energy dependence of the bulk orbital g- factor mo/mc{E) calculated at the 
energy of the IP electron level (dash-dotted line) and its value at the bottom of the conduction band (short dash-dotted line). 



